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We consider the hamiltonian operator associated with planar sec- 
tions of infinitely long cylindrical solenoids and with a homogeneous 
magnetic field in their interior. First, in the Sobolev space T-L 2 , we 
characterize all generalized boundary conditions on the solenoid bor- 
der compatible with quantum mechanics, i.e., the boundary conditions 
so that the corresponding hamiltonian operators are self-adjoint. Then 
we study and compare the scattering of the most usual boundary con- 
ditions, that is, Dirichlet, Neumann and Robin. 
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1 Introduction 

Although the Aharonov-Bohm (AB) effect is a fundamental question in 
quantum physics, and despite the original work on the AB effect has been 
published 50 years ago [H [18] , it is still a very active area of research with 
many open mathematical questions. Here we address some of these ques- 
tions and, our first aim is to try to characterize, in the two-dimensional 
space, all boundary conditions on the (cylindrical) solenoid border S that 
are compatible with quantum mechanics, and whose domains are subspaces 
of the natural Sobolev space % 2 (iS'), where S' is the exterior region of the 
solenoid. 

Our (standard) cylindrical solenoid S has radius a > 0, is infinitely long 
and centered at the origin (with axis coinciding with the z direction), it 
carries a stationary electric current so that there is a homogeneous magnetic 
field B = (0, 0, B) confined to the solenoid interior 5°, and vanishing in the 
exterior region S' . A spinless charged particle of mass m = 1/2 lives in S' 
and has no contact with the magnetic field B. If A is a vector potential that 
generates such magnetic field, that is, B = V x A, the initial (quantum) AB 
hamiltonian operator for such charged particle is given by (with K = 1, c 
and q stand for the speed of light and particle electric charge, respectively) 

#=(p-^A) 2 , p = -iV, domH = Cg°(5'). (1) 

Ahead, we will fix a specific choice of the vector potential. 

This operator H is not self-adjoint and so does not correspond to a phys- 
ical observable; the possible self-adjoint extensions characterize all possible 
physical interaction of the particle with the solenoid border (sometimes ob- 
tained through non-trivial limit procedures) . It is important to note that the 
elements ifi of the domain Cg°(<S') do not touch S (in the sense that ip = 
in a neighborhood of the solenoid) and Hip has a very simple action. Such 
collection of "exotic" extensions might, for instance, include the description 
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of limit situations where singular perturbations (that could also depend on 
time) are slowly turned off. 

It is usually assumed that the domain of the "physical" self-adjoint ex- 
tensions is Dirichlet (i.e., tp = on S), and some theoretical arguments have 
appeared to justify such choice (see [9] and references therein). However, 
here we take an open-minded position and ask about other possibilities of 
boundary conditions. In a general sense, these other conditions correspond 
to the requirement of vanishing of the probability current at the solenoid 
border, and they may model different sorts of interactions between the par- 
ticle and the solenoid. In Section [2] we characterize all of such boundary 
conditions whose domain of the corresponding self-adjoint hamiltonians are 
composed of functions with square integrable first and second derivatives 
(these are rather natural technical conditions in quantum mechanics). The 
operator H above has deficiency indices equal to infinity, and its self-adjoint 
extensions should be compared with the case of solenoid of zero radius dis- 
cussed, for instance, in [21 [7], where the deficiency indices are equal to 2. 

In Section [3] the two-dimensional scattering will be discussed in this 
context, but restricted to the traditional boundary conditions: Dirichlet, 
Neumann and Robin. In fact, the case of Dirichlet was investigated in [20] 
and part of our results are based on the techniques discussed therein. The 
scattering in case of solenoids with zero radius is discussed in [H [21 [21 [IT] . 

For solenoids of radius greater than zero, and the above mentioned tra- 
ditional extensions, we will show that the wave operators exist and are com- 
plete; we also find expressions for the scattering operators and their asymp- 
totic behaviours for low and high energies. Finally, we will find explicitly 
the respective differential scattering cross sections (an important ingredi- 
ent in experiments) and some figures will compare their values. From some 
point of view, the discrepancy among such figures could, in principle, be 
useful for an experimental selection of the boundary condition occurring in 
each situation; in principle it is not obvious which boundary conditions are 
naturally realized in laboratories, and we have found that given two of such 
self-adjoint extensions, it is always possible to find a range of energy so that 
the corresponding scattering cross sections can be distinguished. 

2 Self-adjoint extensions 

In this section we find and characterize an important class of self-adjoint 
extensions of the initial hamiltonian H (see equation (HD), with vector po- 
tential A given, in polar coordinates (r, 8), by A = (A r ,Ag), A r = and 

Aq = , r > a, and $ is the total magnetic flux through the solenoid. 

As will be discussed in Section [31 this hermitian operator has deficiency in- 
dices [8] n+(H) = n^(H) = +oo, and so it has infinitely many self-adjoint 
extensions. 
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The first physical and mathematical point to be addressed is to find 
some self-adjoint operators that may potentially describe the Aharonov- 
Bohm hamiltonian of a charged particle moving in the exterior region S'. 
This particle can not penetrate the solenoid but interacts with S, and the 
boundary conditions that give rise to self-adjoint realizations are the possible 
conditions, from the point of view of quantum mechanics, that may describe 
such interaction with different types of interface materials and limit proce- 
dures. 

We note that in order to classify all such extensions it is necessary to 
make use of some Sobolev spaces T~L S (S) with s < [3]. Furthermore, another 
difficulty is that the domain of the adjoint operator 



is not contained in the space H 2 (S') HE E2 E9, which has proved to be 
natural in quantum problems. Below we shall restrict our arguments to the 
extensions whose domains are contained in H 2 (S'), which will permit us to 
use boundary triples to find these extensions in a quite simpler way (see 
Remark [2J . 

Due to the symmetry of the problem, we shall consider a planar cross 
section; another fact is that the solenoid border in M 3 is not a compact set, 
and so it is not clear how to define the trace operators in the spatial case 
(and we want to avoid this technical point). 

The reader interested only in the final results may go straightly to The- 
orem [2] and the examples that follow this theorem. 

2.1 Boundary triples 

Our way to find self-adjoint extensions of H is via boundary triples (h, pi, P2), 
as described in [8], Chapter 7, so we present a brief account of this technique. 

Definition 1. Let T be a hermitian operator in a Hilbert space T~L. The 
boundary form of T is the sesquilinear mapping T = Tj-* ■ domT* x 
domT* —7- C given by 



Proposition 1. T(^,ry) = 0, for all £, r\ G domT*, if, and only if, T* is 
self-adjoint, that is, if, and only if, T is essentially self-adjoint. 

Boundary forms can be used to determine self-adjoint extensions of T 
by noting that such extensions are restrictions of T* to certain domains T> 
such that r(£,r/) = 0, for all £, 77 G V. By von Neumann theory [8], each 
self-adjoint extension of T is in a one-to-one correspondence with unitary 
operators U : K_(T) — > K + (T) between the deficiency subspaces K± of T; 



domiT = {if, G L 2 (5') : Hijj G L 2 (5')} 



(2) 



:= (T%n) - (d,T*n) , £,7? G domT*. 



(3) 
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denote by T u the corresponding self-adjoint extension whose domains is (T 
denotes the closure of the operator T) 

domT^ = {r, = C + V- ~ Urj- : ( G domT,r?_ G K_(T)}. (4) 

Note that the boundary form T restricted to domT u vanishes. 
Now we recall the concept of boundary triples. 

Definition 2. Let T be a hermitian operator with deficiency indices n_(T) = 
n + (T). A boundary triple (h, pi, P2) for T is composed of a Hilbert space h 
and two linear mappings p\ , p% : dom T* — > h with dense images and so that 

bT T .(S, V ) = ( P i(0,Pi(v)) - (P2(0,P2(ri)) , V^r? G domT*, (5) 

for some constant / i 6 C. (-,■) denotes the inner product in h and 
dimh = n+(T). 

Again, self-adjoint extensions of T are restrictions of T* to certain do- 
mains T> so that r(^, rj) = 0, for all ^,t/ 6 P, and given a boundary triple 
for T, such domains V are related to unitary operators U : h — > h so that 
U P i(0 = P2 (0 and 

<Pl(0,Plfa)> = <P2(0,P2(»/)) = {Upii&UftiT,)) , V£,T, G P. (6) 

The main results we need here are summarized in the following theorem. 

Theorem 1. Lei T be a hermitian operator with equal deficiency indices. 
If (h, pi, P2) is a boundary triple for T, then the self-adjoint extensions ofT 
are given by 

domT u = U G domT* : p 2 (£) = U Pl (0} , 
/'i Vi- ^GdomT^ 7 , 

/or eac/i unitary operator U : h — > h. 

2.2 Boundary triples for the AB operator 

Some self-adjoint extensions of the initial AB operator H will be found. 
It will combine the cylindrical symmetry with the topological property of 
multiply connectedness, that is, the plane with a circular hole, without men- 
tioning the important ingredient of a magnetic potential A with div A = 
in S'. The method can be adapted to other regions with boundaries so that 
the trace construction applies (e.g., smooth and compact boundaries). 

Although a ip(r, 6) G T-L l {S') is not necessarily continuous, it is possible to 
give a meaning to the restriction ip(a,9) = ip\s(9) £ L 2 (£) via the so-called 
trace (more properly, it should be called Sobolev trace) of ip; see ahead. It 



5 



turns out that there is a continuous linear mapping 7 : Cq(IR 2 ) C ^{S') — > 
L 2 (S), 7(<yf(r, 0)) = (p(a,9), that is, there is C > so that 

II7^IIl 2 (S) = \\<p(a,e)\\u(S)<C\\<p\\ H HS'), ^Cj(l 2 ). (8) 

Note that for tp G Cq(M 2 ) the boundary values (p(a,6) are well defined for 
any angle 8. By density, this mapping has a unique continuous extension 
70 : 'H 1 {S') — > L 2 (5), called the trace mapping (see chapters 1 and 2 of [To] 
and also [3|6]), and one defines ip(a,9) := ( 7 qVO(0) f° r an ^ £ 1~L l (S'). 
Similarly it is defined the trace mapping 



7iV>, (9) 



where h* is the normalized vector normal to S pointing to inside the solenoid. 
We shall also make use of the Green's formulae 



Aip(x,y)ip(x,y)dxdy + / Vtfj(x,y)V(p{x,y) dxdy = / yiip -y Q <p da, 

JS' Js 

(10) 

which holds for all i/),tp G T-L 2 {S'), and 

J , fa&y^&'y) ]dxd V+ J / ^{ x ,y)-Q^i x ,y)dxdy = J 7 o^7oV 7o(n-e^)(icr, 



(11) 

— (x,y)(p(x,y)dxdy+ J i/j(x,y)-^-(x,y) dxdy = J 70^ 7oP 7o(n-e^) da, 

(12) 

which hold for all ip,tp G 'H 1 {S'), where tier is the "surface" measure in 
5 (recall that here S is the circle centered at the origin and radius a > 
0), (x,y) G M? and e£ 3 are the unit vectors along the axes x and y, 
respectively. 

Note that the kernel of the trace operator 70 is the Hilbert space 

HliS') := {iP G U\S') : (70^) (0) = ^{a, 9) = 0} , (13) 

which can also be defined as the closure of C§°(5') in H 1 ^'). 

Now we introduce a boundary form V for the initial AB operator (UJ), 
and companion mappings p\ and P2 as well. The boundary form of H , for 
t/},(p G dom //"*, is 

TU>,<p):={H*1> } <p)-W,H*<p), (14) 

and by restricting to those self-adjoint extensions whose domains are con- 
tained in H 2 (S'), Sobolev traces can be invoked. Since divA = 0, the 
boundary form of H is found to be given by 



ry>,(p) = (H*4,<p)-(ii,,H*<p) 

= / {lo4> 71 9 s - 7i"0 7o9?) da - 2i / (70^ 7o(A • n) 70^) da, 
Js Js 
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for all ^,tp£ % 2 {S'). 

By passing to polar coordinates (r,9), the above boundary form T may 
be rewritten as 



2tt 



-dip , . dip 



T(ip,(p) = a I [ ij)(a,0) — (a,6) - — (a,9)cp(a, 



o 



- 2i ip(a, 9) (A • f) (a, 6)<p(a, 9) dO, (16) 



for all tp, (p G H CS'), with and -7— (a, 0) denoting the traces 70X and 

ar 

71X, respectively, for all x G 1~L 2 {S'). 

Now we introduce the boundary triple (h, p\,p2), with h = L 2 (5), acting 
in ^ 2 (S') by Pj : ^ 2 (<S') -> L 2 (S), j = 1, 2, 



/nU-.i r(f/.tf) + / { ^(a,9)-i(A-f)(a,eW(a,e) ) . 
/'2('.') r(r/^)- / :(^(a ) e)-»(A.f)(a,e)V(a,e) j 



(17) 



After a short calculation it follows that 

(2i/a)T(i/),ip) = (pi&),Pi(<p)) h 2 {s) - (P2W,p2(tp)) L 2 (s) , (18) 

for all i/>,(p e % 2 (S')- 

Since A • r = 0, the expressions of pi and P2 are reduced to 



Pl (^) = V>(M) + ^M), 

P2 (^)=V(a^)-^(a^) J 



(19) 



and the vector potential no longer appears in these expressions; note that 
this was possible only due to the cylindrical symmetry of the problem. 

Finally, by applying Theorem [TJ the above constructions permit us to 
conclude the main result of this section: 

Theorem 2. All self-adjoint extensions H u of H, acting in H 2 (S'), are 
characterized by unitary operators U : L 2 (5) — > L 2 (5) so that p2(ip) = 
U pi(ip), that is, 

rU _ L, r- <w2/ C >\ /- rrs.,.,. m , rr^, 



domH u = |V en 2 (S') : (1 - EfyKM) = i(l + U)-^(a,9)j , 
H u ip = H*4>, ip£domH u . 
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2.3 Some self-adjoint extensions of H 

For sake of completeness, in what follows we present some particular choices 
of unitary operators U and the corresponding self-adjoint extensions [8] of 
the initial AB operator H. 

Example 1. If U = — 1, then 

dom^ = G U 2 {S') : z/>(a,0) = 0} = n 2 (S') n Ul{S'), 
H u ^ = H*i} 1 ip£domH u . 

This is the so-called Dirichlet self-adjoint realization, which is usually as- 
sumed to be the physical relevant in the literature |20^ [9] . 

Example 2. If U = 1, then 

domH u = {ip £ U 2 (S') : d^/dr(a,6) = 0}, 
H U 7p = H*ip, ip£domH u . 

This is the so-called Neumann self-adjoint realization. 



(22) 



Example 3. Assume that (1 + U) is invertible. 

In this case, to each self-adjoint operator A : dom^4 C L 2 (5) — > L 2 (5) 
corresponds a self-adjoint extension H A . In fact, first pick a unitary oper- 
ator Ua so that A = - Ua)(1 + Ua^ 1 , domA = rng (1 + Ua) and 
rng A = rng (1 — Ua)', remind of Cayley transform. Now, domH A is the set 
of tp G 7i 2 (S') with "dip/dr(a, •) = Aip(a, •)," understood in the sense that 

(l-U A )ip(a,e)=i(l + U A )^ : (a,e), (23) 

in order to avoid domain questions. Of course the quotation marks may be 
removed in case the operator A is bounded. 

Similarly, for each self-adjoint B acting in L 2 (5) there corresponds a 
unitary Ub, and if (1 — Ub) is invertible, then it corresponds the self-adjoint 
extension H B of H with dom.H B being the set of ip G T-L 2 {S') so that 
u ip(a, •) = Bj^(a, •)," in the sense that 

(l-U B )^(a,e) = i(l + U B )^(a,e). (24) 

The quotation marks may be removed in case the operator B is bounded. 

Note that Example ahead is, in fact, particular cases of this example 
in which A = A4f and B = A4 g are multiplication operators. 

Example 4. U is a multiplication operator. 
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Given a real-valued (measurable) function u{9) defined on 5, put U = 
M. e in(8). If the set {9 : exp(m(0)) = —1} has measure zero, then the function 



1 _ e iu(0) 

is (measurable) well defined and real valued. The domain of the correspond- 
ing self-adjoint extension H u of H is 

domH u = {ip£ U 2 (S') : d^/dr(a, 9) = f(9)ip(a, 9)} . (26) 

Similarly, if {9 : exp(iu(9)) = 1} has measure zero, 

is real valued and the domain of the subsequent self-adjoint extension H 
oiH is 

domH u = {ip G U 2 {S') : ^{a,9) = g(9)d^/dr(a,9)} . (28) 

Special cases are given by constant functions /, g, the so-called Robin self- 
adjoint realization. We shall discuss the scattering for this extension in 
Section [3l 

Remark 1. Since the deficiency indices of H are infinite, there is a plethora 
of self-adjoint extensions of H in the multiply connected domain S' . Some 
of them can be quite unusual and hard to understand from the physical and 
mathematical points of view. 

Remark 2. By using a continuous extension of the trace maps to the dual 
Sobolev spaces H~^ 2 (S) and n-^ 2 (S), in fH one finds references and 
comments to her previous works on all self-adjoint extensions of the laplacian 
in terms of self-adjoint operators from closed subspaces of T-i~ 1 ^ 2 {S) . It is 
possible to follow those works and apply the same technique to find all self- 
adjoint extensions of the initial AB operator H , but we will not describe 
them here since the characterizations are rather abstract, involve spaces not 
usual in quantum mechanics, and they require a length construction that 
is not so clean as the extensions we have found in Theorem [D ( which also 
includes the traditional self-adjoint extensions we are most interested in). 



3 Scattering 

In this section we study and compare the scattering for the Robin self-adjoint 
realizations (which includes Neumann and Dirichlet as particular cases) of 
the initial AB operator H. 
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3.1 Scattering theory: a brief account 



Now we briefly recall results from scattering theory, based mainly on |20|, [5] , 
focused on what we want to do in the next sections. For details and a 
thorough study of the mathematics of scattering theory see [W\ |2"T] . 

Consider a system in (non-relativistic) quantum mechanics whose states 
£ are unit vectors in a Hilbert space 7~L and whose time evolution is gener- 
ated by a self-adjoint operator h acting in H, and let Hq denote the free 
hamiltonian acting in T-Lq. The question is whether the states e~ lht ^ are 
scattering states, i.e., if there are free states £± G Hq so that 



£ - e iht Je- iHot Z± 



(29) 



vanishes as t — > ±oo. A comparison mapping J : Ho — > H, which is a 
unitary operator (or just a bounded one), is sometimes conveniently intro- 
duced. 

Definition 3. The wave operators are the strong limits 



W± := s- lim e M Je~ lHot , (30) 

t— i>±oo 



if they exist. 



Recall that the wave operators W± are said to be complete if rng W± = 
^p(/i) _L , where H p (h) denotes the closure of the subspace spanned by the 
eigenvectors of h. 

The vectors £± and £ satisfy the relation £ = W±£± and the wave op- 
erators W± : domW± — >• rng W± are partial isometries. Thus, W±Wijl 
are orthogonal projections onto rng W±, and restricted to these subspaces 
W± = W± . Furthermore, 

C+ = S^, (31) 

where S := W^W_ is the so-called scattering operator or S-matrix. 

The physical system we consider is a scattering of particles off a cylin- 
drical obstacle in the plane of points x = (x,y), and the Hilbert space is 
H = L 2 (5'). Moreover, we are also interested in the scattering away from 
a short range continuous potential V(x,y), which we also assume that it 
is spherically symmetric, that is, V(x, y) = V(r), r = \{x,y)\; short range 
means that there are constants C, R > so that 

\V(x,y)\<^, \fr>R, (32) 



for some 5 > 0. Let 



-A + V(r) and H = p\ + p|, (33) 
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acting in the position space H and momentum Hq = L 2 (M 2 ), respectively, 
where p = (pi,P2), P 2 = \p\ 2 and the comparison operator J is the inverse 
Fourier transform T . 

In the time-independent scattering theory one solves the time-independ- 
ent Schrodinger equation hip = p 2 ip for the incoming tp_(x,p) and outgoing 
<p + (x,p) functions, which are reduced to a plane wave (ft(x,p) = e lx ' p for 
\x\ — > oo (these are solutions of the Schrodinger equation of the free particle). 
One has the following connection 

(W±C)(x) = ^ J v±(x,p)C(p)dp, C e H . (34) 

By employing polar coordinates, both in space position (r, 9) and in the 
space of momenta (k, 9'), one obtains 

h= -w-^-hw +v(r) and Ho = k2 > (35) 

respectively. One then considers the asymptotic behaviour of the solution 
(incoming wave function) of the time-independent Schrodinger equation 

ip-{r,6;k,6') ~ e ikrcos( - e -^ + f(k,6 - 9')^-, r ^ oo, (36) 

where / is the scattering amplitude, and so the differential scattering cross 
section is 

^y k ,6) = \f(k,9)\ 2 . (37) 

Physically, this quantity measures the probability density of an incident 
particle, after interaction with the target (scatterer center), i.e., a scattered 
particle to be found within of a cone around (k,9). 

Now we recall how to find the asymptotic behaviour of ip^ and the 
scattering amplitude. If , J n denotes the Bessel function of first kind of 
order n, one has 

oo oo 

e ikrcose = J- i^J lml (kr)e ime , f(k,9)= ]T f m (k)e im8 , (38) 

m=— oo m=— oo 

and by (f36j) . the asymptotic behaviour of J n (r) [16J, 



and recalling that cost? = (e^ + e )/2, one obtains 

oo 

p_(r,0;A;,O)~ £ 



/ ^in / 4—ikr / g— i7r/4 j 

^nkr) 1 / 2 + W27rA;r) 1 /2 + ^7 



(40) 
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On the other hand, if - solves the Schrodinger equation and is regular at 
the origin. By using separation of variables, 



tp-(r,0;k,0) = a m (k)^ m (r,ky me , (41) 

m=— oo 

where (f m is a solution for radial Schrodinger equation of angular momen- 
tum m, 

d 2 1 d m 2 Tr . A l9 , a s 

"55 + IT + V(0) «. = *>,>, (42) 

which is regular at the origin. Since V is short range, it follows that the 
behaviour of ip m for r — > oo is given by 

ip m (r, k) ~ ( -^-] cos ( kr - hm\ir - y + 5 m (k)) , (43) 



% 7rfcr/ \ 2 4 

and herein the phase shift 6 m was introduced. The phase shift is a measure of 
the argument difference to the asymptotic behaviour of the solution Ji m i(kr) 
to the radial free equation that is regular at the origin. 

The comparison of the asymptotic behaviour of ip- given by (|40p . with 
the one given by (|4ip and (|43p . gives us the important relations 

a m (k) =i\™\ e i& ^ k \ (44) 

g 2i(5 m (fc) _ i 

fmik) = 7^)^' (45) 

and 

1 oo 

/(M) = ^ E (e 2 ^«-l)e^, (46) 



m=— oo 



which formally expresses the scattering amplitude in terms of the phase 
shifts 6 m . Later on, in our applications to the extensions of the AB hamilto- 
nian, these relations will be considered from the point of view of distribution 
theory. 

Now a bit of the time-dependent approach again. Decompose the Hilbert 
spaces H and Ho in subspaces h m , ho, m > with corresponding projections P m , 
Po,m, respectively. For example, 

(P mV )(r,0) = —-^r, m (r), (47) 

where 

rim(r) = j^j- 2 fj e-* me V (r, 6)d6 (48) 
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are the components of rj with angular momentum m, belonging to the Hilbert 
space % r = Ly dr ([0, oo)). The subspaces h m are invariant under the hamil- 
tonian h and the subsequent restriction is given by 

d 2 1 d m 2 T r , , , J „. 

h - = -^--r^ + ^ + V ^ (49) 

in T-L r . For the inverse Fourier transform T one has 

oo 

(T()(r,6) = (2tt)- 1/2 e ime (Fm( m )(r), (50) 

m=— oo 

with T m : Tlk H r denoting the unitary operator 

roo 

(T m iP)(r) = i\ m \ J\ m Akr)rl>{k)kdk, (51) 



where Tik = L| dS .([0, oo)). Correspondingly, we consider a sequence of wave 
operators from Tik into H r 

W±, m = s- hm e^^e"^. (52) 

t— >±oo 

Since W±m exist and are complete in each sector m, one has the following 
relation with the time-independent approach 

/•oo 

(W±, m V)M = * H / ip m (k,r)e^ 5 ^i;(k)kdk, (53) 



so that the corresponding S'-matrix in the sector m, S m : - > Hk, 

S m = W; >m W_, m , (54) 
is given, after some calculations, by 

{S m m) = { e ^ m{k)m ' \ = k y , (55) 

that is, /Sm, is simply the multiplication operator by e 2%&m ^ on Ti^. 
The wave operator in Definition [3] can now be written as 

oo 

(W±C)M) = (2^r 1/2 £ e^(W±, m C m )(r). (56) 

m=— oo 

Similarly, the scattering operator S 1 : %o - ^ %o is found to satisfy 

"oo r2ir r2ir 



/'OO />Z7T />Z7T 

(C,^) = (C,0 + (2vr)- 1 E / / / ^^aM) 

m =-oo J V J JO 

x f(jfe,0') ( e 2l< M fc ) _ ^ dOdO'kdk, 



(57) 



13 



and by equation ([35]) , we conclude that 



f'OO pIty p'2ti 

(CS0 = (CO+ / / C(M)£(M') 

Jo Jo Jo , co s 

( ik \ l ' 2 (58) 

x - f(k,e -e')d6de' kdk, 

for all C,£ G where (•, •) denotes the inner product in Ho, and so 

"(5-l)(M)=(^J /(M)," (59) 

understood in the sense described above, which is the correct relation be- 
tween / and S. 

3.2 Robin self-adjoint extensions 

In this subsection we describe the self-adjoint extensions of the initial AB 
operator ([T]) for which we will study scattering. We choose some of the self- 
adjoint extensions that preserve angular momentum since we are considering 
that the subspaces h m are invariants under H. In addition, the extensions 
described below are the most common and studied in the literature when 
borders are considered. Note that there are extensions that do not preserve 
angular momentum. 

In order to simplify expressions, we will take c = q = 1 in our initial 
hermitian operator ([1]), 

H = (-iV - A) 2 , (60) 

which acts in a subspace of the Hilbert space T~L = L 2 (S'), and recall that 
the vector potential A, in polar coordinates (r, 8), is given by A = (A r , Aq), 

with A r = and Aq = - — , r > a. This operator can be written in polar 
coordinates as 

H--— --— — (%— -aV (61) 
dr 2 r dr r 2 \ 80 J ' 

where a = — 3>/(27r). Without loss of generality, consider < a < 1, and 
a = means that no magnetic field is present. 

The next step is to construct the self-adjoint extensions of H we are 
interested in. By making the polar decomposition % = H^He, where %® = 
L 2 dr [a,oo) and He = L 2 [0,27r], we obtain a sequence of formal restriction 
operators to h m 

d 2 1 d jm + a) 2 

ttm+a — --T2 ~T~ H 2 ^° l > 

dr z r dr r z 
in To achieve our goal we need to turn these operators into self-adjoint 
ones acting in They are not essentially self-adjoint on C X) (a, oo) for any 
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m + a. To see this, note that the potential term (m + a) 2 /r 2 is bounded in 
W", hence we need only to consider the differential operator 

h - ^ 1 d + 1 (63) 
dr 2 r dr 4r 2 

But using the unitary operator U : Hf. — > L^a, oo), given by (Uip)(r) = 
r 1 / 2 ^(r), the operator h\ becomes 

h 2 = UhU- 1 = - Jl (64) 

which is not essentially self-adjoint on C§°(a, oo) since the functions u±{r) = 
e -e ±Z7T / 4 r g L^ r [a, oo) and satisfy h^uziziu = 0; in other words, its deficiency 
indices are re_(ii2) = n+(ii2) = 1. Since this holds for all m, it justifies the 
assertion that n±(H) = oo in Remark [TJ 

However, we can find all self-adjoint extensions of h.2 [8], which are well 
known and given by 



dom 



h\ = fa G U 2 [a, oo) : tp(a) = A^'(a)} , h\^ = h 2 ij), (65) 



for each A G M U {oo}. 

Thus, we have the corresponding self-adjoint extensions of H^^.^ 

dom Hi +a = U- 1 (dom h^, H^+J, = H m+a ^, (66) 

that is, 



domi^ +Q 



fa = r ~ 1/2 ip : tp G H 2 [a, oo) and if; (a) = A?// (a)} 

u G U- 1 (n 2 [a,oo)) : r^uML^ = X^-W^uir) 
v dr 



(67) 



Therefore, the boundary conditions that characterize the Robin self-adjoint 
extensions of H m+a are given by (2a — A)u(a) = 2a\u'(a). If A ^ 2a, then 

domil^+a = \ uGU- 1 (n 2 [a, oo)) : u(a) = ^-W (a) \ . (68) 
\ 2a — A J 

We shall denote these Robin self-adjoint extensions by that is 

dom#4 +Q = {u G C/- 1 (^ 2 [a,oc)) : «(a) = An'(a)} , 

A ^ 

where A = 2aX/(2a — A). 
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Note that upon integrating by parts it follows that if A > 0, then 
(H^+aU, u) > 0, for all u G dom.H^ a+a , that is, the self-adjoint opera- 
tor H^ l+a is non-negative and, therefore, cr(ff^ +a ) C [0,oo); from now on 
we assume that A > 0. 

Since each sector is invariant under H, the Robin self-adjoint extension 
of full operator H is given by 

H X = H x m+a ® 1, (70) 

and note that this extension is a special case of Example U and the choices 
of constant functions /, g guarantee that such self-adjoint extensions pre- 
serve angular momentum. The principal cases are for A = and A = oo, 
which correspond to the well-known self-adjoint extensions of Dirichlet and 
Neumann, respectively. 



3.3 Scattering for the Robin extensions 

In this subsection we study the scattering for Robin self-adjoint realizations 
of the initial AB operator H introduced in equation (pQ). We find the scat- 
tering operator S, we prove the existence of wave operators and that they 
are complete, and we also obtain explicit expressions for them. In addition, 
we determined the scattering amplitude and hence the differential scattering 
cross section for such extensions. 

We underline that there is no magnetic field in case a = 0, and the 
scattering is sole due to the presence of the solenoid; in fact, the Aharonov- 
Bohm effect is noticed by comparing the results for with this reference 

case a = 0. 



3.3.1 Scattering operator 



Assume initially that the wave operators exist and are complete; under such 
conditions, we shall write out the scattering operator. The solution to 



d 2 Id (m + a) \ 



dip 



(71) 

with the linear combination ld — A— vanishing at r = a, where A = -, 

dr 6 2a - A 

is given by 

<Pm(k, r) = G^{k, a) (A r | m+Q | (ka) - \NL +a , (ha)) J\ m+a \ {kr) 



(ka)) N\ 

m+a\ 

(kr) , (72) 
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with G^ n (k,a) to be determined by imposing condition (|43p . By using the 
asymptotic behaviour of J v and N u for r — > oo, we obtain 



irkr 



1/2 



sin I fer |m + a| 7r 



(j| m+Q |(A;o) - AJ|' m+Q | (fca) 



(73) 



By comparing the above expression with (|43p we have 



Gm(&, a) 



cos fer m + a 7T 

2 1 1 4 



1 



sin ( At - ^\m + a\ir -j) (j\ m+a \{ka) - A J| m+a | (ka 



cos [ kr — - \m\ it — — + 5^(k, a) 



(74) 



that is, 



7T 



cos I kr — - \m + a\ tt — — + 6\ \ = cos [kr — - \m\ tt — — + 5 m (k, a 



7T 



with 0^ so that 

Q N lm+al (ka) - XN; m+al (ka) 
cos 6\ = — , sin d\ 



(75) 

J\rn+a\(ka) - A JL +Q;| (fco) 



(76) 



with 

and, therefore, (|i"3j) is satisfied if 



^/ (N\ m+a \(ka) - \N> m+al (ka)y + (j| m+Q |(/ca) - AJ^ +a| (fca))' 



(77) 



(78) 



Note that D never vanishes. In fact, suppose that D = 0. Then N\ m+a \ (ka) — 
AA ? "| / m+Q |(A;a) = and J| m+Q |(A;a) — A «7j' m+a i (ka) = 0. So, it follows that 
J\m+a\(ka)N( m+a \(ka) - N\ m+a \(ka)J' m+a ^(ka) = 0. But this is a contradic- 
tion with the Wronskian W r=a [Ji m+Q ,i (kr), N\ m+a \ (kr)] = 2/(ira) ^ 0. 
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Now, comparing the arguments of the cosines above, it is found that the 
phase shift 5^(k, a) is given by 



S^(k,a) = A m (a)+9 X , 



(79) 



7T 



\m + a\). Therefore, the scattering operator 5^ m : 
7ik — > Hk for the Robin self-adjoint extension is 



where A m (a) = '- (\m\ 



= 2i<5* (fc,a) = _ e 2iA m (a) 



{J\ m +a\(ka) - iN\ m+a \(ka)) - A (j( m+a i(ka) - iN',(ka) 



(J\ m +a\(ka) + iN\ m+a \(ka)) - A f J[ m+a \{ka) + iN'Aka) 



(80) 



that is, 



a 2iA m (o) 



h\ 2) , .(ka) - XH\\ ,(ka 

\m+a\ y I \m+a\ v 



(2)' 



,(ka) - \Hl ±r , ,(fco) 



(1)' 



|m+o| 



where H^'^ 2 \x) = J v (x) ± iN u (x) are the Hankel functions [T| I10j. 



(81) 



Remark 3. Note that for A = we get the Dirichlet case and recover the 
expression of the scattering operator found in [201; and, if we choose A = oo 
we obtain the scattering operator for the Neumann case, namely, 



S: 



= 2iA m (a) l m +"l V ' 



(82) 



3.3.2 Asymptotic behaviours of the scattering operator 

Now we describe the asymptotic behaviour of the scattering operator for 
different self-adjoint extensions for both ka — > oo and ka — > 0, and also 
compare the results. This is done from the asymptotic behaviour of Bessel 
functions. 

We begin with the behaviour for ka — > oo. For this we recall that 

(\m ~\~ c\\ \ 
— <^|m+a| + l(^°)^ H ~ J\m+a\ 

(ka) , (83) 

and thus, its behaviour for ka — > oo is given by 
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and so 

2 \ 1/2 



J'm + a\(ka) ~ - (^) sin (ka - \m + a\\ - |) k 

1 /2 

|m + al / 2 \ ' / 7T 7r\ 

H — — cos [ka — \m + a\ . 

a \irka V 2 4/ 



(85) 



Similarly, 



/ _|_ Q< 

N \m+a\( ka ) = ( -N\ m+a \ + i{ka)k + N\ m+a \{kc 



(86) 



and its behaviour for ka — > oo is given by 



o \ !/ 2 

2 \ /, , , vr 7T 



cos ( A;a — |m + al ) k 



irka J V 2 4 

|m + a I / 2 \ ^ 2 / 7r 7r 

H — — sin [ka - \m + a — - — 

a V vrfca / V 2 4 



(87) 



Now, by considering A 7^ 0, the asymptotic behaviour of Bessel functions 
and their derivatives given above, we get for ka — > 00 

<-iA / -I ™ — 2ika+iir/2 

(/ca) 2 A 2 + 2iA(A|m + a| — a)/ca - A 2 |m + a| 2 + 2A|m + a|a - a 2 



(fca) 2 A 2 + A 2 |m + a| 2 - 2A|m + a\ + a 2 



(88) 



and, since the term in square brackets is approximately 1 in this case, we 
have 

Sa, m ~ (- l) m e - 2ifea +^/ 2 , (89) 

for fca — >■ oo. This expression coincides with the Neumann case, i.e., the 
scattering operator for the Robin case acts as the Neumann case for large 
energy, independently of A, provided that A ^ 0. However, it differs from 
the Dirichlet case 

qV i -i\m —2ika—i-n 12 ((\(W 

for ka — > oo. 

Summing up, for very large energies the scattering operator does not 
distinguish different Robin extensions (i.e., < A < oo) from the Neu- 
mann case A = oo, but it has a different behaviour from the Dirichlet 
case A = 0. In order to try to understand such behaviour intuitively, let 
us informally consider the perhaps simplest situation, that is, the "free" 
unidimensional reflection from a barrier at the origin x = with wavefunc- 
tion ip(x) = A(k)s'm(kx) + B(k)cos(kx), at least near the origin; Dirich- 
let and Neumann boundary conditions impose that B = and ^4 = 0, 
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respectively, whereas Robin condition ^(0) = Aip'(Q) imposes the energy 
relation \Xk\ 2 = \B(k)\ 2 /\A(k)\ 2 , and for large energies k — > oo one has 
|-A(fc)| "C and in this region the system behaviour becomes similar to 

the one dictated by the Neumann condition. 

On the other hand, taking into account that the behaviour of the Bessel 
functions for ka — > [16] are given by 



and 



J\ m+a \(ka) ~ [(l/2)A;a]l m+a l/r(|m + a| + 1), 
at ii \ T(\m + a\) 



(91) 



(92) 



^^[(l/2)ka]\ m + a \ , 

we get the following behaviour for the scattering operator m for ka — > 



5 A ~ 



d? - d 2 (ka/2) A \ m+a \ n 2d 1 d s (ka/2) 2 \ m+a \ 



+ i 



d 2 + d 2 (fca/ 2 ) 4|m+a| 



. Q d 2 -d 2 {ka/2f\ m+a \ , 



d? + d 2 (W2) 4|m+Q| 

2did 3 (fca/2) 2 l m+Q l 



(if + d 2 (fca/ 2 ) 4|m+a| d 2 + d 2 (fca/2) 4 l m + Q l 

with (3 = 7r(|m| — \m + a|) and the coefficients 



(93) 



di 
d 2 
d 3 



-r(|m + a|)/vr - A[2T(|m + a\ + 1) - |m+ a\T(\m + a\)]/(ira), 
r(|m + a\ + 1)" 2 (1 + A 2 |m + a| 2 /a 2 - 2A|m + a\/a), 
T(\m + a\ + 1) _1 (1 - A|m + a\/a), 



(94) 



are independent of k. For m = and a = we have 



"Sao ~ 



41n(fca)' 



+ A 



2 _ 7T 2 (fca/2) 2 
aln(fca) aln(fca) 2 



\2 ( l-f{ka/2Y 
A I a 2 ln(fca) 2 



tt 2 , x / 2 , ^ 2 (W2) 2 \ , X2 / 1+^ 2 (W2) 4 \ 
41n(fca) 2 A ^aln(fea) aln(fca) 2 J ^ A \ a 2 ln(fca) 2 ^ 



1 + 



+ i 



7T 



1 + A 



a ln(fca) 



+ 



2(fca/2) 2 



+ A 



2 2(fca/2) 2 
a 2 ln(fca) 



ln(fca) 



1 + 



41n(fca) 2 



+ 



A \ a\n(ka) ~*~ a\n(ka) 2 J A \ a 2 ln(fca 



/2)4 

r 



(95) 



for fca —7- 0. 

We observed that, for very small energies fca — >• 0, the sole scattering 
operator is not able to distinguish the Robin, Dirichlet and Neumann self- 
adjoint extensions of the initial AB hamiltonian ([I]). Furthermore, this 
occurs both for the case with field (a ^ 0) and for the reference case without 
field (a = 0). These behaviours were numerically recovered. 
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3.3.3 Wave operators 

Now, we prove that in fact the wave operators for the Robin self-adjoint 
extensions exist and are complete. In addition, we obtain an explicit ex- 
pression for them. We begin with a lemma that will be used in the proof of 
Theorem [3l 

Lemma 1. Let A and B be self-adjoint operators and U a bounded operator 
so that AU = UB, assuming that the compositions are well defined. Then 

e~ iAt U = Ue- iBt . (96) 

Proof. For each £ S domB so that U£ £ dom^, let u(t) = Ue~ lBt £ s and 
v(t) = e~ iAt U£, for all tel. We want to show that u(t) = v(t), for all 
t € R. If w(t) = \\v(t) - u(t)\\ 2 , then 

II = Jt^ ~ V{t) ~ U{t)) = ~ Jt [v{t) ~ U{t)] 

= 2Re [(-t) {(e- iAt Ui,Ae- lAt UC) + (Ue~ iBt £, AUe~ iBt £) 
- 2Re{Ae- lAt U£,Ue- lBt £))} =0 

(97) 

because what is in square brackets is real since A is self-adjoint. Therefore 
w(t) is constant. Since w(0) = 0, it follows that w(t) = 0, for all t G R. □ 

Let P a : U r -> be given by (P a ip)(r) = (x[o,oq)V')('")> that is ) is 
the orthogonal projection operator onto 

Theorem 3. The wave operators 

Wi am = s- lim e^-'Pfl^e^' (98) 

t— »±oo 

exis£ and are surjective isometries from Hk onto T~lf.. Explicitly, for ip E Hk, 
they are given by the expressions 

( a m A (r) = fH lim [* ^(k,r)e Ti5 ^ k ' a ^{k)kdk. (99) 

Proof. We consider only the case a m \ the proof for W+ Q m is similar. 
Since the proof is rather long, we divide it in three steps: 
I s * Step: Define the candidate for the limit operator, and show two equali- 
ties. 

2 nd Step: Show that the wave operator a m exists and satisfies ([99|) . 
3 rd Step: Show that the wave operator a m is a surjective isometry. 
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1 Step: "Define the candidate for the limit operator, and show two 



equalities." Let us define an operator Uh a m by 



ip(k)kdk, 



(100) 



with ip eHk and supp ip C (e, i?). This operator is well defined by Holder's 
inequality. By using that 

J H {kr) = (2/tt) 1 / 2 cos (kr - \v\tt/2 - it / A) / {kr) 1 / 2 + 0((A;r)- 3 / 2 ) (101) 

and 

N\ u \{kr) = (2/vr) 1 / 2 sin (kr - \v\k/2 - ir/A)/(kr)^ 2 + 0((fcr)~ 3 / 2 ) , (102) 



it is found that Uh a m ^> G and it satisfies 



17* ^ Wo <^A(i2,e)W Wfc , 



(103) 



where K\(R,e) depends only on R and e. 

Now we check that Ub. a m il) is actually in the domain of H^ l+a and 



(104) 



In fact, let u G domi7^ +a . Then, by Fubini's theorem, we can write 



/ (^+a«)(0(^-,a,mV')(0^r (105) 
J a 



"f[( 



_d 2 __lj_ + (m + a) ^ _ (r) , H 



dr 2 r dr 



x / ^(k,r)e l5 ™ {k > a ^{k)kdk 



rdr (106) 



■H 



1 f 6 



<5^(fc,a) 



^(fc) 



/•oo j 

(m + a) 2 —^(k^rya^rdr 

J a r 



J 

■J a 



^m( fc - r ) 1^2 + ^ 



fedfc, (107) 



and integrating by parts the second term in square brackets, we obtain 

-R 



(f 2 Id (m + a) 2 \ \ 
d^~rd? + 7 2 ' if m (k,r)u(r)rdr 



kdk; (108) 
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d * 

note that there are no boundary terms since (p^kja) = A j^ (k, a) and 
u(a) = Xu'(a). Therefore, 



dr 



R 



and, again by Fubini, we find that 



fc 2 ^(A;, r)u(r)rdr 



kdk, 
(109) 



R 



u(r)i H I ip^(k,r)e i5 ^ k ' a) k 2 i){k)kdkrdr, 

(110) 



that is, 



(111) 



where (•,•) denotes the inner product in Then U_ am ^> S domH x l+a 
and 

^ +a ^, a ,J = U^ m k 2 ^. (112) 
Finally, apply Lemma Q] to conclude 



e -i H ik+ a t U x ib = U X e~ ik2t ib 



(113) 



2 S'iep: "Show that the wave operator exists and satisfies ([99]) ." As- 
sume that ib G Co°(e,i?). By using the conclusion of the first step (i.e., the 
last equality above), one can write, 



e iH ^P a T m e- ik \ - U x ib 



P a T m e- lk \ - U^ a<m e- ik2t iP 
(P a T m - U^Je-* k \) (r 



rdr 



(114) 
(115) 



(116) 



J ]m] (kr)e- ik2t i/j(k)kdk 



R 



i |m| / dkkip^(k,r)e 



ijj(k)kdk 



rdr (117) 



■\m\ 



R 



-ik 2 t„ 



[J\m\(kr) - <p m (k,r 



A (h r \J S m( k >°') 



kdk 



rdr, (118) 
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and after some calculations with the asymptotic behaviour of the two func- 
tions in brackets above, we obtain 



R 



-ik 2 t 



Akr 



K?(ka) 



(fcr)V2 



{ka)0 ((fcr)- 3 / 2 ) 



kdk 



rdr, 
(119) 



with and K% are functions of class C°°. Finally, using the inequality 



+ 9\r <2 / r + 2Nr, we get 



< 2 



R 



e~ ik V(fc)i^i (ka) 



Skr 



(kr) 1 / 2 



kdk 



rdr 



(120) 



+ 2 



ik \{k)K^(ka)0 ((kr)^ 2 ) kdk 



rdr. 



We discuss each term on the right side of the last inequality separately. 
By replacing e~ tk t+tkr with (— 2ikt + ir)~ Y dke~ lk t+lkr \ integrating by parts 
and using the dominated convergence theorem to estimate the first term, it 
is found that it vanishes as t — >■ — oo. 

For the second term, let h(kr) = O ((/cr) -3 / 2 ) , then h(kr) = M{kr) 
x(kr)~ 3 / 2 , with M(kr) a bounded function. So, by Riemann-Lebesgue 
lemma and dominated convergence theorem, 



R 



-ik z t 



%b(k)K^(ka)0 ((kry 3/2 ) kdk 



rdr 



R 



-ik z t 



ip{k)K$(ka)M(kr)k- 1/2 dk 



r- 2 dr -> 0, (121) 



as t — > — oo. Then 



0, 



(122) 



as t — > — oo, and since Co°(0,oo) is dense in T~Lk, it follows that the wave 
operator a m exists and satisfies (f99j) . 

3 rd Step: "Show that the wave operator a m is a surjective isometry." 
To show that the wave operator is an isometry we take £rlk with compact 
support and check 



W il) 

' ' -,a,mr 



lim 

t — y — oo 



P a T m e- lk \ 



lim 

t— >— oo 



(123) 
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To prove that rng W_ 



it suffices to show that its adjoint is an 



isometry, because the kernel {0} = N ((W* am )*) = (rng W^ iQjm ) , and 
so rng Q m = Since am (W^ a m )* is the orthogonal projection 
onto rng Q m , which is closed, and since 



i' m ' lim 



R-^oo 
,-R 



R 



y x m {k,r)e iS ^ ((W* Q ,J"V) (k)kdk 



R 



V^{k,r)e i5 ™ {k ' a) { (-i)H / ^(M)e~*" (fc,a ty(s)a<fe ) fcdfc 



lim 







G^(k,a)D^(ka,kr) / G^ n (k,a)Dl(ka,ks)ib{s)sdskdk, 



(124) 



and recalling that G^(k, a) = — , with 



I? = ^ (JV„(*;a) - XNl(ka)) 2 + (J v (ka) - XJ^ka)) 2 , (125) 

and = |m + a|, one has 

£>*(fca,i/) := [#„(fca) - AiV^(fca)] J„(y) - [j„(fca) - A.#(fca)] A^(y), (126) 

and, in order to conclude the theorem, it is enough to prove the following 
lemma. 

Lemma 2. Let ib G C§°(a, oo) and z/ > 0. T/ien 



ip(r) 



lim 



1/R 



D 2 



D^(ka,kr) / D^(ka, ks)ip(s)sds kdk. (127) 



In fact, this lemma implies that 



(W*JV (*) (r)=^(r), 



(128) 



for all ^ G Co°(a, oo), and since this set is dense in it follows that 



rng W_ 



and the theorem is proved. 



□ 



In the following we present the proof of Lemma [2j 

Proof. Since the Wronskian of J v (z) and N v {z) is equal to 2/(ttz) [IB], that 
is, W Z [J U ,N U ] = 2/(vrz), one has W r [J u (kr), N u (kr)] = 2/(vrr). Now, we 
consider the boundary value problem 

(E - H u )ld = ib, a < r < oo, lim .El > 0, 

(\ \>(\\ (129) 
(p[a) — Ac/? [a) = (J. 
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The Green's function 



ui(r)u 2 (s) 



g(r\s) = < 

k W s [ui,u 2 

is the solution to the auxiliary problem 



W s [u 1 ,u 2 
ui(s)u 2 (r 



a < r < s, 
s < r < oo, 



(130) 



where 



[E — H u )g = 5(r — s), a < r < oo, 
g(a) - Xg'(a) = 0, 

ni (r) = \N u {E l / 2 a) - XN' v {E l l 2 a)\ J v (E^ 2 r) 
- JuiE^a) - A4(£ 1/2 a)l N u {E l / 2 r) 
= D x v (E l l 2 a,E l ' 2 r), 



(131) 



(132) 



is the solution to (E — H u )u = that satisfies the boundary condition at 
r = a, and 

u 2 (r) = H^(E 1 / 2 r), (133) 

is the solution to (E — H u )u = that satisfies the boundary condition at oo, 
and the superscripts (1) and (2) correspond to ImE > (with ImV^E > 0) 
and ImE < (with Im^/E < 0), respectively; W s [iii,U2] is the wronskian 
of the solutions u\ and u 2 at the point r = s, and in this case one has 



W s [ui,u 2 ] = — (u 2 (a) - A^(a) ) . 
its \ dr 



(134) 



Write Re ■= (E — H v ) 1 for the resolvent of H v at "energy" E, so that 
the solution (i?£f/0( r ) to problem (|129|) is given by 



(R E ^)(r) 



g(r\s)tp(s)ds 
ui(s)u 2 (r) 



a W s [ui,u 2 ] 



ui(r)u 2 (s) 



ip(s)ds 



7T 

2 



u 2 (a) - A^(a) 
dr 



-i 



W s [ui,u 2 ] 

H (lH2) {E l/2 r) I D ^ E l/2 a ^ E l/2 s) 



xij(s)sds + D x {E 1 / 2 a,E 1 / 2 r) / H®'® {E 1 l 2 s)^{s)s ds . 

J r 

(135) 

Recall now the Stone formula [8] for the spectral projection of H v onto 
the interval [a, b], 

1 f b 

X[a,b] {Hu) = s- lim — I (R x _ i5 - R x+i5 )dx. (136) 

5->0+ ZTTl J a 
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If one writes E- = x — i5 and E + = x + i5 for the "energy" E with Im.E < 
and ImE > 0, respectively, then 

[X[l/R,R)( H M ( r ) = , lim + f l&E-Mr) - (R E+ ^)(r)] dx. (137) 

Now, substitute the above expressions for the resolvent operators and use 
the dominated convergence theorem (take into account that the functions 
J u and N v are continuous and bounded), after some manipulations and 
simplifications we obtain the expression 

[x[i/R,R](H„)ip] (r) 



__ 1 f R Dl{x l ! 2 a,x l / 2 r) 

" 2 Ji /R (N u (x l / 2 a) - \Nl(x l l 2 a)f + {J u (x l / 2 a) - XJ^x^a)) 2 

POO 

x / D*(x 1/2 a,x 1/2 s)ip(s)sdsdx. (138) 

J a 

Finally, use the change of variable x 1//2 = k, so that dx/2 = kdk, to get 

/R ^ roo 
j^D^(ka,kr) / D^(ka,ks)i)(s)s ds k dk, 

(139) 

which is a.e. equal to the function defined by the integral on the right side 
of (|127|) . On the other hand, X[a,b]{Hu) = for [a, b] C (— oo,0) since H v is 
a positive operator (see Theorem 8.3.13 in [8]), and so cr(H v ) C [0, oo) and 
H u has no eigenvalues. □ 



3.3.4 Scattering amplitude and cross section 

In this subsection we calculate the scattering amplitude and differential scat- 
tering cross section for the Robin self-adjoint extensions, and some compar- 
isons will be made in the next subsection. 

However, first we recall what was done in [20] to determine the scattering 
amplitude f a for the case of a solenoid of radius zero, and with Dirichlet 
condition at the origin; since we will make use of such results. In the case 
of radius zero, in each sector of angular momentum m, one has 

A m (a) = ^(\m\-\m + a\), (140) 

for the phase shift, which is a function only of a, and so the corresponding 
scattering operator is 

^(a) = ( ej, rn>~a 
e , m < —a 
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Then, the Fourier coefficients of f a in the expression (|46p has constant mod- 
ulus and do not vanish as \m\ — > oo; so the scattering amplitude f a is seen 
as a distribution. To obtain the correct expression of the amplitude f a , note 
that the scattering operator S a on Hq is an integral operator and, by using 
the above expressions, in [20J it was found that 

(S a £)(k, 9) = s a (9 - e')de', (142) 







with 



AO) - W ) oos( xa ! + i 8 J^plpv f-J-^ ) . ill:;. 



where PV denotes the principal value (recall that here < a < 1). These ex- 

1 /2 

pressions and the relation (S — l)(k, 9) = [jj^) f(k, 9) imply the following 
expression for the scattering amplitude 



/2tt\ 1/2 



, . „, .sin(7ra) __, / 1 
<5(0)[cos(7ra) - 1] + i — i — -PV 



ir \e ie - 1 

Now, if 9 the distribution f a is represented by the function [20 



(144) 



sin(Tra) e^/ 2 
^ 6 >~ (2^)1/2^(0/2)' (M5) 

and so the differential scattering cross section in this case is 

do~\ /, „■> 1 sin 2 (7ra) 

which agree with the expressions found by Aharonov and Bohm [3] and also 
by other authors, for example in [T3|. Thus we will continue looking at the 
scattering amplitude as a distribution, which will be conveniently calculated 
from the Fourier series. 

Remark 4. We observe that in \1$ it is advocated that there should be no 
5(9) in the above expression (|144p for the scattering amplitude, and that f a 
should be restricted to (|145p . - this causes a controversy with references [Lflj 
and WU§ . In any event, since we do not consider the forward direction 9 = 
in our comparisons of the scattering due to different self-adjoint extensions 
(i.e., our main goal in the next section), we are able to keep away from such 
controversy. 

Now we turn to our Robin extensions and positive radius. For a = 0, 
that is, no magnetic field, the expression (|46p gives for scattering amplitude 
associated with H x 

^ oo 

^ / m = — fYi 
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and since 

we obtain 

/o(M) 



„2i**(fc,0) 



Hr\(ka) - XHr\ (ka 



\m\ 



H^(ka)-XH^(ka) 



irik 



1/2 



E 



m=—oo 



J\m\(ka) - XJ'(ka) 



irnO 



H^(ka)-XH^(ka) 



(148) 



(149) 



Note that for fixed ka and k ^ 0, the series above is convergent since its 
coefficients are fast decaying as \m\ — > oo, due to the well-known behaviour 
of Bessel functions. Thus, in this case /q is represented by a function and 
therefore the differential cross section is given by 



2 

irk 



m=— oo 



(150) 



On the other hand, again by (|4*6|) . the scattering amplitude associated 
with iJ^, with non-zero magnetic field, that is, < a < 1, is given by 



and since 



(2irik) 



3 2iA m («) 



_. oo 



(151) 



m=— oo 



XH, 



(2)' 



(Jfeo) 



#i L i(ka)-XHl\ .(ka 



(152) 



we obtain 



(27ri£;) 1/2 



E 



m=— oo 



3 2iA m («) 



(2) 

|m+a| 



(Jfea) 



^|SU|(*») 



(i) 



1 e 



imd 



(153) 



Now, let n G Z be fixed and change variable m' = m + n in the summation 
index. Then m = m' — n and since A m (a) = (7r/2)(|m| — |m + a|), we obtain 



-in 9 



fa(k,6) 



{2mk) 1 l 2 



E 



m '= — oo 



= 2i<5 m , (a-n) 



(-ir 



flf^ .(Jfeo) - XH\ 2) ' ,(ka) 

\m'+a— n\ y ' \m'+a— n| v ' 

,(ka) - XhP',. .(jfeo) 



1 e l 



(154) 
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which can be written as 

ti(k,9) = (-l)« e - i ^_„(* ! ,e)+(27r/ifc) 1 / 2 [(-l) n -l]<y(e), n G Z. (155) 

Thus, the differential cross section for the Robin self-adjoint extension of 
the initial AB hamiltonian is given by {9 7^ 0) 



c£cr\ ' 



1 



E 



m=— oo 



= 2iA m (a) 



JET, 



(2) ( 
m+a| ' 



za)-XH^ (ka 



,(ka) - XH>\ ,(k 

\m+a\ v ' \m+ot\ v 



(1)' 



\m+c 



+ 1 e 



iraS 



, (156) 



which is periodic in a with period 1. This is a justification for the restriction 
< a < 1. It is convenient to write 



ti(k,0) = f a (k,9) + f r ,x(k,e), 



(157) 



where / a is the scattering amplitude of the case of radius zero a = with 
Dirichlet condition at the origin, which was discussed above, 



f a (k,0) = {2mk)~ l l 2 Yj ^ ; " A;! ' ,! 

m=— a 

and with f r \ given by 



1 e 



irnO 



(158) 



fr,\{k, 9) 



■nik 



1/2 



E 



m=— oo 



„2iA m (a) 



a) J\ m+a \(ka) - XJ' lm+al (ka) 



imd 



Xka)-XH^ ] ' ,(ka) 

\m+a\ \ ' |m+a|\ / 



(159) 

By the same argument presented above, the series for f r \ is convergent, and 
/«(M) is given by CHS]). 

Therefore, the differential cross section for the Robin extension with 
parameter A, for k 7^ and 9 ^ 0, is given by 



sin(7ra) e 



-i0/2 



(2vriA:) 1 /2 s in(0/2) 



3 2iA m (a)_ 



imd 



B\ X \ :(ka)-XH, 



r(l)' 



(ifeo) 



(160) 



Again, A = corresponds to the Dirichlet case, whereas A = 00 to the 
Neumann boundary condition. 
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3.4 Scattering comparison 



In this section we present some figures and comments to illustrate and com- 
pare the scattering results obtained in the previous subsections. In the 
figures, we have fixed the value of the solenoid radius to a = 1. Due to the 
symmetry of the differential cross sections as function of 6, the correspond- 
ing plots are presented only for < 6 < it. Recall that the scattering in 
case a = is simply due to the solenoid of non-zero radius, and one notices 
the Aharonov-Bohm effect by comparing this case with the scattering for 
different values of a (in particular for non- integer a). 

In the following, we collect the expressions for the scattering operators 
for Dirichlet, Neumann and Robin extensions, respectively, 

v _ cos/3 [N\ m+a \{ka) 2 - J\ m+a \(ka) 2 ] - 2 sin (3 J\ m+a \(ka)N\ m+a \(ka) 
a ' m ~ N\ m+a \(ka)* + J\ m+a \(ka) 2 



+ i 



sin/3 [iV| m+a |(£;a) 2 - J\ m +a\(ka) 2 ] + 2 cos /3 J\ m+a \ (ka)N\ m+a \ (ka) 



N\ m+a \(ka) 2 + J\ m+a \{ka)< 



(161) 



cos (3 



' N \ m+a \(ka) 2 ~ J[ m+ 4 ka A ~ 2 s ™P J \ m+a \(k*)N' lm+al (ka) 

l +tt iH 2 + l +a ,H 2 



\m+a\ 

sin/3 I TV' M ,(ka) 2 - J! M Akaf 



+ i 



\m+a\ ' 

2 cos B J( , ,(ka)N! , ,(ka) 



^ m+al (ka) 2 + J' lm+al (ka) 2 



(162) 



and 



3-A 

J a,m 



COS (3 



(^V|m+a| (ka) - AA f | / m+Q | (ka)} - {j\ m+a \ (ka) - \J[ m+ct \ (kc 



N\ m +a\ (ka) - AAf | ' m+Q| (ka)} + (j\ m+a \ (ka) - A J^ +q| (ka)} 
2 sin/3 (j\ m+a \(ka) - \J[ m+a \(ka)} (N\ m+a \(ka) - AJVj^ (fca)) 
( N \m+a\ (ka) - AiV' m+a | (ka)} + (j| m+a | (fco) - A J^ +q| (fca)) 
sin /3 (iV| m+Q | (fca) - AA^ | ' m+a| (fca)) - ( J| m+Q | (ka) - A J[ m+a \ (ka)} 



(N\ m+a \ (ka) - \N' {m+al (ka)} + (j\ m+a \ (ka) - \J[ m+a \ (h 
2cos/3 (j\ m+a \(ka) - \J[ m+a \(ka)} (N {m+al (ka) - AiV | ' m+a| (ka 



(^|m+a| (ka) - A^' m+Q| (fca)) + (j| m+Q | (fca) - AJ' m+a | (ka)} 



(163) 
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and recall that /3 = 7r(|m| — \m + a|). 

Figure [Q presents the real parts of scattering operators for the above 
three extensions, in a case with a / 0. Note that, for high energies, the 
curve of the Robin scattering operator approaches the curve of the Neu- 
mann case, and it is evident the phase difference between the Dirichlet and 
Neumann cases; this agrees with the theoretical results of Subsection 13.3.21 
For very low energies (ka — > 0), this figure illustrates what we have said in 
the last paragraph of Subsection 13.3.21 about the behaviour of the scattering 
operator (since cos/3 = 0), that is, for low energies the scattering operator 
is very similar in all cases we have considered. Similar results hold when no 
magnetic field is present, i.e., a = 0. 

For each of such self-adjoint extensions, we have numerically checked 
(with plots) that the scattering operators, with non-zero magnetic fields 
(i.e., for any a / 0), approach the corresponding scattering operators with 
no magnetic field (i.e., a = 0) for high energies (not shown). Hence, given 
one of those self-adjoint extensions, for high energies the scattering operator 
is not able to discern the presence of magnetic field inside the solenoid or 
not. We note that such behaviours of the scattering operators were found 
to be independent of the values of m, < a < 1, and A > 0. 

Now, we consider the important concept of differential cross section, 
in the case of cylindrical solenoids of positive radius a > 0, and for the 
Dirichlet, Neumann and Robin extensions. The respective expressions we 
have obtained are, for 9 ^ 0, 





2 



(164) 





2 



(165) 



and 





2 



(166) 
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For high energies, we have found that the differential cross section of 
Neumann and Robin cases are very close to Dirichlet for each given < a < 
1, except in a neighborhood of 9 = and 9 = 2tt. Figure [2] shows those 
curves for a = 1/2. 

Figure [3] shows the differential cross section of the three extensions in 
terms of the "energy" k, for the case with non-zero field, represented by 
a = 1/2, fixed angle 9 = ir/2 and A = 1; note the different behaviours for 
high and low energies. 

For k — > 0, in the case with field (a ^ 0) and positive radius, we have 
found that the differential cross sections for the three cases have the same 
behaviour, which is approximately given by the differential cross section of 
the case with zero radius (|146p and Dirichlet condition at the origin. See 
Figure HI 

For intermediate energies the differential cross sections of the extensions 
differ significantly, as illustrated in Figure this seems interesting, since it 
is an explicitly distinction among different boundary conditions. 

Finally, we mention that for small A (for example, A = 1/10), the differ- 
ential cross section for the Robin extension approaches the values obtained 
for the Dirichlet case, and when we choose A large (for example, A = 10) the 
values for the Neumann extension are virtually recovered. This is certainly 
expected. 

4 Conclusions 

With respect to the mathematical problems related to the traditional mag- 
netic AB setting, that is, the one associated with an infinitely long solenoid, 
in this work we have based our investigations on two cornerstones. First, we 
have considered the more realistic case of a solenoid of positive radius a > 0; 
and second, we did not take for granted that the boundary conditions on 
the solenoid border S is Dirichlet (although there are physical insight [20J 
and mathematical arguments that support this choice |9j). 

The boundary conditions that are physically compatible with quantum 
mechanics are those that define self-adjoint extensions of the initial AB 
hamiltonian (TTJ). We have characterized all such self-adjoint extensions 
whose domains are contained in the natural Sobolev space 7-L 2 {S'); this was 
done via boundary triples, and our main contribution was the inclusion of 
the vector potential in the operator action, by taking into account the sym- 
metry of the problem, and a gauge choice as well, to simplify expressions. 

The important cases of Dirichlet, Neumann and Robin are among the 
self-adjoint extensions we have characterized via boundary triples, and the 
next step was to study the scattering for these self-adjoint hamiltonians; 
such study was based on [20], where the particular case of Dirichlet boundary 
condition was considered. For some parameter ranges, that is, < A < oo, 
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we have proven that the wave operators are well defined and complete; 
furthermore the hamiltonian is positive and has no eigenvalues. We remark 
that for negative values of A one can not discard the presence of eigenvalues 
(see, for instance, Exercise 7.3.3 in [8]), and so bounded states could emerge 
from the Robin boundary condition; this is an interesting possibility we 
think it is worth investigating. 

Then we have explicitly calculated the scattering operators and subse- 
quent scattering cross sections, and they were our natural physical quantities 
used to compare different self-adjoint extensions. Note that the scattering 
cross section is a distribution in general, but for the scattering angle 9 ^ 
it is represented by a continuous function (k ^ 0). 

For high energies, we have found that the scattering operator for the 
Robin case is similar to the Neumann one, but different from the Dirichlet 
case. On the other hand, for low energies the behaviour of the scattering 
operator is independent of these self-adjoint extensions. Such results hold 
for each fixed < a < 1. 

With respect to the differential cross section, for "intermediate energies" 
its behaviour depends significantly on the choice among the three self-adjoint 
extensions we have considered. 

To finish, we underline that, in general, our scattering results depend 
on the magnetic field parameter a, and this is actually a confirmation of 
the presence of the AB effect (when < a < 1) in different self-adjoint 
extensions! 
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Figure 2: Differential cross section as function of 6 in the case with field 
(a = 1/2), with a = 1, k = 30 and A = 1/10. 
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Figure 3: Differential cross section as function of k in the case with field 
(a = 1/2), with a = 1, 9 = n/2 and A = 1. 
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Figure 4: Differential cross section as function of 6 in the case with field 
(a = 1/2), with a = l,k = 1/10 and A = 1. 
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Figure 5: Differential cross section as function of 6 in the case with field 
(a = 1/2), with a = 1, k = 3/2 and A = 1. 
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